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(57) ABSTRACT 

This invention is a ground flutter testing system without a 
wind tunnel, called Dry Wind Tunnel (DWT) System. The 
DWT system consists of a Ground Vibration Test (GVT) 
hardware system, a multiple input multiple output (MIMO) 
force controller software, and a real-time unsteady aerody- 
namic force generation software, that is developed from an 
aerodynamic reduced order model (ROM). The ground flutter 
test using the DWT System operates on a real structural 
model, therefore no scaled-down structural model, which is 
required by the conventional wind tunnel flutter test, is 
involved. Furthermore, the impact of the structural nonlin- 
earities on the aeroelastic stability can be included automati- 
cally. Moreover, the aeroservoelastic characteristics of the 
aircraft can be easily measured by simply including the flight 
control system in-the-loop. In addition, the unsteady aerody- 
namics generated computationally is interference-free from 
the wind tunnel walls. Finally, the DWT System can be con- 
veniently and inexpensively carried out as a post GVT test 
with the same hardware, only with some possible rearrange- 
ment of the shakers and the inclusion of additional sensors. 
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DRY WIND TUNNEL SYSTEM 

STATEMENT OF GOVERNMENT 
INVOLVEMENT 

5 

This application involves an invention made with United 
States Government support under a Small business Technol- 
ogy Transfer (STTR) Program, Phase I entitled, “GVT-Based 
Ground Flutter Test without Wind Tunnel,” having Contract 
Number: NNX08CD34P awarded by NASA Dry den Flight 10 
Research Center (National Aeronautics and Space Adminis- 
tration’s Dry den Flight Research Center, Edwards, Calif. 
93523). 

FIELD OF THE INVENTION 1 5 

The present invention relates generally to a method by 
which software tools and hardware systems simulates flight 
tests/wind tunnel tests on an aircraft, part of an aircraft, or any 2 o 
other flexible structure to search the aeroelastic/ aero ser- 
voelastic (AE/ASE) instabilities of the test structure. More 
particularly, the present invention relates to a method which 
uses software and hardware components for the simulation of 
a flexible structure’s AE/ASE characteristics at various flight 25 
conditions. 

BACKGROUND OF THE INVENTION 

Prior to flight testing an aircraft structure or any part of a 30 
flexible structure, flutter wind-tunnel tests are required to 
search for AE/ASE instabilities of the structure’s configura- 
tions. Flutter wind-tunnel testing is an expensive and compli- 
cated process because it involves the design and fabrication of 
a scaled-down aeroelastic structural model, which is a scaled- 35 
down version of a real structure, and also because wind- 
tunnel time is costly. Flutter wind-tunnel tests are necessary 
because if a flutter instability is not found prior to the flight 
test then major instabilities may occur during the flight tests 
which can make the structure extremely unstable, hence 40 
endangering the aircraft structure and passengers. Unfortu- 
nately, creating a scaled-down version of a real structure may 
introduce discrepancies in the structural characteristics of the 
scaled-down structure when compared to the real structure. 

For example, the control surface actuators with correspond- 45 
ing stiffness and damping, are difficult to represent accurately 
in the scaled-down structure. Also, the incorporation of the 
accurate modal damping in the scaled-down structural model 
is almost impossible. These discrepancies between the actual 
structure and the scaled-down versions can lead to uncertain- 50 
ties in the measured aeroelastic instability boundary. In addi- 
tion, during flutter wind-tunnel tests the wind-tunnel walls 
interfere with the test results. 

It is very costly to include structural nonlinearities such as 
friction, free-play, etc. in the scaled-down structural model. 55 
And even when these structural nonlinearities are included in 
the scaled-down version of the structure they will never be 
identical to the real structure. These structural nonlinearities 
can have a major impact on the aeroelastic characteristics of 
the aircraft. The measurement of the coupling between the 60 
flight control system and the aeroelastic system to search for 
ASE instabilities of the aircraft is another critical design 
requirement that can be accomplished by performing a flutter 
wind-tunnel test in the presence of flight control laws. How- 
ever, performing this aeroservoelastic (ASE) measurement in 65 
the wind tunnel is also a very expensive process, thus the ASE 
analysis largely relies on flight tests. 
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If the flutter wind-tunnel test was able to use the real- 
structure then the control laws would automatically be 
included, since the real structure has control laws included in 
it to be able to fly. With the use of a real -structure there would 
be no scaling discrepancies between the tested structure and 
the actual structure. The impact of the structural nonlineari- 
ties on aeroelastic stability can automatically be included. If 
the flutter test could be conducted without a wind-tunnel then 
the unsteady aerodynamics generated computationally would 
be interference- free from the wind-tunnel walls. 

In order to be able to conduct a flutter test without a wind- 
tunnel, a ROM (Reduced Order Model) of the unsteady aero- 
dynamic model is needed. This unsteady aerodynamic ROM 
represents an aerodynamic transfer function that inputs the 
physical structural deformation and outputs the aerodynamic 
forces. Thus, this requirement immediately rules out the 
Computational Fluid Dynamics (CFD)-based ROM’s 
because all CFD-based ROM’s involve some type of modal 
approach that assumes the structural mode shapes are known. 

SUMMARY OF THE INVENTION 

The object of this invention is to provide a method for 
accurate calculations of an aircraft structures’ aeroelastic 
instabilities at a comparatively low testing outlay. 

The term “Dry Wind Tunnel System” (“DWT System”) 
refers to a system that consists of a Ground Vibration Test 
(“GVT”) hardware system, a data acquisition system, a real- 
time unsteady aerodynamic force generation software mod- 
ule and a multiple input multiple output (MIMO) force con- 
troller software module. 

As used herein, the term “Ground Vibration Test hardware 
system” (“GVT System” or “GVT Hardware” or “GVT 
Hardware System”) refers to a test that includes state-of-the- 
art equipment including devices such as aircraft soft- support, 
structural excitation, vibration sensors, and data acquisition 
and analysis tools. The GVT Hardware used by the DWT 
System is essentially a ground vibration test system which 
consists of vibration sensors, shakers, and load cells. The 
vibration sensors are devises such as accelerometers, linear 
variable differential transformer’s (LVDT), and linear veloc- 
ity transducer (LVT). The accelerometers measure the accel- 
eration of the structure at each sensor location. The LVDT’s 
measure the displacement of the structure at each sensor 
location. The LVT’s measure the velocity of the structure at 
each sensor location. The shakers are a structural excitation 
device, in which they exert forces onto the structure at each 
shaker location. Each shaker can exert forces different from 
the other shakers. The load cells are devices which measure/ 
monitor the amount of force being exerted by the shakers onto 
the structure. 

The real-time unsteady aerodynamic force generation soft- 
ware is a reduced order model (ROM) of the unsteady aero- 
dynamics that has a real-time computational efficiency for 
unsteady aerodynamic force generation at different flight 
conditions. This computer software program reads the struc- 
tural responses as its inputs, and the desired unsteady aero- 
dynamic forces at each shaker’s location are computed by 
passing the structural response input to the reduced order 
model of the unsteady aerodynamics. 

The ROM that the real-time unsteady aerodynamic force 
generation software uses is created by the GVTOPT Module 
of the ZAERO Software System (ZAERO is ZONA Technol- 
ogy, Inc.’s commercial software product for aeroelastic/ aero - 
servoelastic analysis). The basic requirement in formulating 
an unsteady aerodynamic ROM for the dry wind tunnel test is 
that it contains only the aerodynamic characteristics of the 
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configuration and does not involve any structural character- 
istics. This is because the coupling between the aerodynamics 
and structures is achieved only during the real-time test; 
thereby the structural characteristics cannot be known prior to 
the test. 

Prior to the DWT System test, the ROM matrices, com- 
puted by the GVTOPT Module at each flight condition, are 
saved in a look-up table database. During the DWT System 
test, the ROM matrices saved in the look-up table can be 
directly called by the real-time unsteady aerodynamic force 
generation software to generate the real-time unsteady aero- 
dynamic force for flutter testing. 

The data acquisition system is such a device that connects 
hardware devices (such as sensors and actuators) and com- 
puters. The data acquisition system collects the analog signals 
measured from sensors, converts it into digital signals, and 
sends the converted digital signals to the computer for analy- 
sis. Visa versa, the data acquisition system can also covert the 
digital signals coming from the computer into analog signals, 
and send them back it to the hardware devices. 

The MIMO force controller software is performed to 
ensure the forces exerted to the flexible aircraft structure by 
the shakers are the desired unsteady aerodynamics forces, 
which are computed by the real time unsteady aerodynamic 
force generation software at each shaker location. Further- 
more, the MIMO force controller is applied to minimize the 
interference among the shakers and reduce the nonlinear 
effects due to the nonlinear characteristics of the shaker 
dynamics. The MIMO force controller software can be 
designed using the standard robust control design technique 
such as TO and \i tools. To facilitate the MIMO force control- 
ler design process, a system identification technology is 
implemented to estimate a MIMO shaker dynamics, instead 
of using the analytical tool for a complicated mathematical 
shaker model. Therefore, an experiment is performed by 
exerting a set of sweep sinusoidal signals to the shakers to 
excite the structure. The input sweep sinusoidal signals and 
the output forces signals on the structure measured by load 
cell sensors are used for model estimation of the shaker 
dynamics. 

The DWT System concept would be particularly useful as 
a pre-flight testing effort to identify any aeroelastic and aero- 
servoelastic instability on new and/or modified aircraft Inher- 
ent structural nonlinearities such as friction and freeplay are 
notoriously difficult to model properly in linearized analyses, 
but would naturally be present in the DWT System testing as 
it is carried out on the actual structure. The ground flutter test 
using the DWT System operates on the real structural model, 
therefore, unlike in the flutter wind-tunnel test no scaled- 
down structural model is involved. 

Also, the impact of the structural nonlinearities on the 
aeroelastic stability can automatically be included. In addi- 
tion, the aeroservoelastic characteristics of the aircraft can be 
easily measured by simply including the flight control system 
in the loop. The unsteady aerodynamics generated computa- 
tionally would have no interference from the wind-tunnel 
walls. The DWT System takes full advantage of the existing 
GVT Hardware setup in that the required DWT System’s 
setup, for flutter testing, uses similar hardware as the GVT 
System of the same aircraft structure, only with some possible 
rearrangement of the shakers and the inclusion of additional 
sensors. Once setup, the DWT System will perform the 
ground flutter test and requires only a link between the soft- 
ware programs and the GVT Hardware. Testing using the 
DWT System would also be useft.il as a post-flight testing 
procedure to resolve any discrepancies between the analysis 
and flight test results. The DWT System concept is applicable 
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to a broad range of test structures, from components such as 
flexible wings to full aircrafts. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

FIG. 1 demonstrates the interaction between the GVT 
Hardware System, the data acquisition system, and the com- 
puter system in which the MIMO force controller and the 
real-time unsteady aerodynamic force generation software 
to are uploaded. 

FIG. 2 illustrates a structure with sensors. 

FIG. 3 illustrates the GVT Hardware System. 

FIG. 4 illustrates the Dry -Wind Tunnel System with all of 
the component and how they interact each other. 

1 5 FIG. 5 illustrates the real-time unsteady aerodynamic force 
generation software process. 

FIG. 6 illustrates another real-time unsteady aerodynamic 
force generation software process. 

20 DETAILED DESCRIPTION OF THE INVENTION 

The present invention may be described herein in terms of 
various functional components and various methods/steps. It 
should be appreciated that such functional components may 
25 be realized by any number of structural components, hard- 
ware, and/or software configurations to perform the specified 
functions. For example, the present invention may employ 
various electronic components, hardware, and software, 
which can be suitably configured for various intended pur- 
30 poses; such as testing a real aircraft structure or part of a real 
aircraft structure to predict the instabilities of said structure. 

FIG. 1 demonstrates the data flow between the GVT Hard- 
ware System 100, the data acquisition system 200, and the 
computer 300, in which the real-time unsteady aerodynamic 
35 force generation software 400, and the MIMO force control- 
ler 500 are installed on. 

FIG. 2 illustrates a structure 101 with sensors 102 on it. The 
sensors 102 are things like accelerometers, which measure 
the oscillation at each sensor’s location; linear velocity trans- 
40 ducer (LVT), which measure the velocity at each sensor’s 
location; and linear variable differential transformer’s 
(LVDT), which measure the structure’s displacement at each 
sensor’s location. 

FIG. 3 illustrates the GVT Hardware System 100. The 
45 GVT Hardware System includes sensors 102a, 102/?, and 
1 02c, placed on a structure 101 . Shakers 103 exert forces onto 
the structure. The load cells 104 measure the amount of force 
the shakers are exerting onto the structure. The sensors 102c, 
102 /?, and 102c, measure the displacement, velocity and 
50 acceleration at each sensor location. 

FIG. 4 illustrates the Dry -Wind Tunnel System process. 
During the flutter test, one of the shakers 103 creates a small 
impulse as the initial disturbance to excite the structure 101. 
The sensors 102c, 102 /?, and 102c are used to measure the 
55 displacement, velocity and acceleration at each designed 
locations of the test structure. The data acquisition system 
200 transfers the data measured by the sensors 102c, 102/? 
and 102c to the computer 300. Within the computer the infor- 
mation from the sensors 102c-102c, goes to the real-time 
60 unsteady aerodynamic force generation software 400. The 
real-time unsteady aerodynamic force generation software 
400 uses the data from the sensor 102c-102c, to calculate the 
desired forces the shakers 103 should exert at each shaker 103 
location onto the structure 101 to search AE/ASE instability 
65 of the structure 101. After the desired forces have been cal- 
culated by the real-time unsteady aerodynamic force genera- 
tion software 400, the GVT Hardware System 100 starts 
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exciting the structure to search for the potential AE/ASE 
instability of the structure 1 01 . The load cells 1 04 measure the 
force exerted by the shakers 103. The data acquisition system 
200 passes the load cell 104 measurements to the computer 
300. Within the computer 300 the measurements are passed to 5 
the MIMO force controller 500. The MIMO force controller 
500 acts to enforce the forces exerted by the shakers 103 onto 
the structure 101 are same as the desired forces generated by 
the real-time unsteady aerodynamic generation software. 
Specifically, the MIMO force controller 500 uses the load cell 10 
104 measurements to monitor if the forces exerted by the 
shakers 103 onto the structure 101 are the same as the desired 
force generated by the real-time unsteady aerodynamic force 
generation software 400. If any discrepancy between the 
measurement of the load cell 104 and the desired force gen- 15 
erated by the real-time unsteady aerodynamic force genera- 
tion software 400 occurs, the MIMO force controller 500 
generates a correct excitation force signal to the shakers 103 
to excite the structure 101. 

FIG. 5 demonstrates the real-time unsteady aerodynamic 20 
force generation software’s process for determining the cor- 
rect force that should be applied by the shakers 103 (FIG. 4) 
onto the structure 101 (FIG. 4) to search for the AE/ASE 
instabilities of the structure 101. The real-time unsteady aero- 
dynamic force generation software is a reduced order model 25 
(ROM) of the unsteady aerodynamics that has a real-time 
computational efficiency for unsteady aerodynamic force 
generation. All the variables presented in FIG. 5 are pre- 
computed via the GVTOPT module of the ZAERO software 
system at different flight conditions, and a lookup table for 30 
these variables can be generated to incorporate variation of 
the flight conditions. 

The basic requirement in formulating an unsteady aerody- 
namic ROM for the dry wind tunnel test is that it contains only 
the aerodynamic characteri sties of the configuration and does 35 
not include any structural characteristics. This is because the 
coupling between aerodynamics and structures is achieved 
only during the real time test; thereby the structural charac- 
teristics cannot be known prior to the test. 

In fact, this unsteady aerodynamic ROM represents an 40 
aerodynamic transfer function that inputs the physical struc- 
tural deformation and outputs the aerodynamic forces. Thus, 
this requirement innnediately rules out the CFD-based ROM 
because all CFD-based ROM’s involve some type of modal 
approach that assumes the structural mode shapes to be 45 
known. On the other hand, the unsteady aerodynamic panel 
methods such as the Doublet Lattice method DLM, ZONA6, 
and ZONA7 do readily generate an aerodynamic influence 
coefficient (AIC) matrix that relates the structural deforma- 
tion (x) to the aerodynamic forces (f a ) which reads: 50 

V a {<*)}=qJAIC(to)]{X} (l) 

where q^ is the dynamic pressure, and w is the oscillating 
frequency indicating that the AIC matrix is formulated in the 
frequency domain. 55 

Closely examining the AIC matrix, one can see that each 
coefficient in the matrix represents the response on a receiv- 
ing panel to a unit input on a sending panel. For a number 
of aerodynamic panels in the aerodynamic model, the size of 
the AIC matrix is N^xNp. Thus, the AIC matrix truly repre- 60 
sents an aerodynamic transfer function that contains only the 
aerodynamic characteristics of the configuration. The AIC 
matrix generated by the linear unsteady aerodynamic meth- 
ods can accurately produce aerodynamic forces at the sub- 
sonic and supersonic Mach numbers. For transonic flows, the 65 
transonic AIC matrix can be generated by using an unsteady 
transonic method called ZTRAN. ZTRAN uses an overset 


field-panel method to solve the linearized transonic small 
disturbance equation with the steady background flow being 
imported from a high fidelity CFD code such as a Navier- 
Stokes solver. The accuracy of the ZONA6 and ZONA7 
methods for subsonic and supersonic Mach numbers and the 
ZTRAN method for transonic Mach number has been dem- 
onstrated on many test cases. 

To obtain a time-domain aerodynamic transfer function, 
requires transforming the frequency-domain AIC matrix into 
time domain. This can be achieved using the minimum state 
method that fits the AIC matrices at a set of frequencies in to 
a rational function approximation which reads: 


ifait)} = q a 


Do] + y~ DlF + y2 D2F 2 |{M + 


(2) 


goo [D] 



mx} 


where s is the Laplace variable, L is the reference length, V M 
is the free-stream velocity. [A 0 ], [AJ, [A 2 ], [D], and [E] are 
the resulting matrices from the minimum state method, and 
[R] is a diagonal matrix with the diagonal coefficients being 
the N z assigned aerodynamic lags. Thus, the size of [R] is 
N z xN z . Usually, no more than ten aerodynamic lags are suf- 
ficient to obtain an accurate rational function approximation. 

The vectors {X} and {f^(t)} in Eq. (2) are the structural 
deformations and aerodynamic forces at the aerodynamic 
panels, respectively. Therefore, for aerodynamic panels, 
the size of the vectors {X} and {at)} is N^x 1 . However, in the 
GVT the structural deformations are measured at the sensor 
locations and the aerodynamic forces are produced at the 
shaker locations. To transfer these structural deformation and 
aerodynamic forces from the aerodynamic panels to the sen- 
sor and shaker locations requires an interpolation procedure 
which can be accomplished by using the spline method such 
as the infinite plate spline method, the thin plate spline 
method, and the beam spline method that jointly generate a 
spline matrices for such displacement and force transferal. 
For N s accelerometers, the displacement spline matrix reads: 

{Ahvpx 1 = [Gs\NpxN S {X S en }a^x 1 (3) 

where {X^ ew } is the structural deformation at the N s sensors, 
and [GJ is a displacement spline matrix that transfers the 
structural deformations from the sensor locations to the con- 
trol points of the aerodynamic panels. 

For shakers, the force spline matrix reads: 

{fshak\N A y. 1 = [Gj\n A x/Vpl/UiVpx 1 (4) 

where f shak is the aerodynamic forces at the shaker locations, 
and [Gy\ is the force spline matrix that transfers the forces at 
the control points of the aerodynamic panels to the shaker 
locations. 

Substituting Eqn.(3) and Eqn.(4) into Eqn.(2) yields a 
reduced-order model of the unsteady aerodynamics shown in 
the following equation: 


[_ L l _ , 

{fshak } = goo | Do] + -p— Dl F + -pj D2F ■Dsen} + 


( 5 ) 


goo [D] 
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where [A 9 ],[A 1 ],[A 2 ]=[G y ][[A 0 ],[A 1 ],[AJ][GJ with size of 
N^xN s , D=[Gy\ {D} with size of N A xN L , E=[E][GJ with size 
of N z xN s -. 

In a typical GVT set-up, the numbers of sensors and shak- 
ers are typically small. 5 

Accordingly, the size of the matrices in Eqn.(5) is small, 
implying that Eqn.(5) can be computed very efficiently to 
generate aerodynamic forces at the shaker locations in real- 
time. 

The real-time unsteady aerodynamic force generation soft- l o 
ware 400 reads the structural responses 601 (FIG. 4) from the 
sensors 102a- 102c (FIG. 4) as its inputs. The inputs 102a- 
102c shown in FIG. 5 correspond to the output of the sensor 
102a-102c in FIG. 4, and which are LVT (displacement), 
LVDT (velocity) and accelerometer sensors. The desired 15 
force 603 computed by the real-time unsteady aerodynamic 
force generation software 400 is passed to the shakers 103 
(FIG. 4) through the data acquisition system 200 (FIG. 4). 

FIG. 6 demonstrates another process of the real-time 
unsteady aerodynamic force generation software 400. In this 20 
process it uses the displacement 102a and acceleration 102c 
as inputs, and outputs the desired force 603 to the shakers 103 
(FIG. 4) through the data acquisition system 200. The veloc- 
ity 102 b is calculated by the integration of the acceleration 
102c. 25 

Both processes demonstrated in FIGS. 5 and 6 can be used 
to determine the desired force. The selection of which process 
to use is based on what information and which sensors are 
available. 

What is claimed is: 

1. A Dry Wind Tunnel (DWT) method that is used to 
predict the flutter boundaries of a structure without using a 
wind tunnel, said DWT method comprising the following 
steps: 

(a) providing: 

(i) at least one shaker, for shaking such structure, 
coupled to such structure at at least one shaker loca- 
tion; 

(ii) at least one computer conununicatively coupled to 
said at least one shaker and storing and operable to 
execute a force controller software for generating 
force control signals for controlling said at least one 
shaker; 

(iii) at least one displacement sensor mounted on such 
structure; 

(iv) at least one velocity sensor mounted on such struc- 
ture; 

(v) at least one acceleration sensor mounted on such 
structure; 

. 50 

(vi) at least one load cell mounted between said at least 
one shaker and such structure; and 

(vii) a data acquisition system communicatively coupled 
to said at least one displacement sensor, said at least 
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one velocity sensor, said at least one acceleration 
sensor, said at least one load cell, and said at least one 
computer; 

(b) loading a real-time unsteady aerodynamic force gen- 
eration software stored in and executable on said at least 
one computer that is able to compute initial real-time 
unsteady aerodynamic force on said structure at said at 
least one shaker location; 

(c) controlling, using said at least one computer, said at 
least one shaker to exert an initial impulse on such struc- 
ture; 

(d) receiving, in said at least one computer, initial data from 
said at least one displacement sensor, said at least one 
velocity sensor, and said at least one acceleration sensor 
via said data acquisition system, responsive to said ini- 
tial impulse; 

(e) computing test inputs for controlling said at least one 
shaker, responsive to said initial data and using said 
real-time unsteady aerodynamic force generation soft- 
ware; 

(f) controlling said at least one shaker, responsive to said 
test inputs and using said force controller software, to 
shake such structure; 

(g) measuring, via said at least one computer: 

i) structural response of said structure using said real- 
time unsteady aerodynamic force generation software 
and structural response data from said at least one 
displacement sensor, said at least one velocity sensor, 
and said at least one acceleration sensor; and 

ii) said force exerted on said structure, using said force 
controller software and said at least one computer, 
responsive to data from said at least one load cell; 

(h) computing, using said at least one computer running 
said force generation software and responsive to said 
measured structural response as input, an unsteady aero- 
dynamic force at said at least one shaker location; 

(i) applying said computed unsteady aerodynamic force to 
said structure through said at least one shaker; 

(j) controlling, using said force controller software, said 
applied forces to be the same as said real-time unsteady 
aerodynamic forces computed using said real-time 
unsteady aerodynamic force generation software; 

(k) monitoring said structural response data in real time; 

(l) if said structural response data indicates a decay motion, 
incrementally changing said inputs to create incre- 
mented test inputs to said real-time unsteady aerody- 
namic force generation software; 

(m) repeating steps (f) to (1) until said structural response 
data indicates a divergent motion; 

(n) designating a condition between said decay motion and 
said divergent motion of such structure as a predicted 
flutter boundary of such structure. 



